Statistics of quantum transport in chaotic cavities with broken time-reversal 

symmetry 
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The statistical properties of quantum transport through a chaotic cavity are encoded in the traces 
Tn = Tr[(tt^)"], where t is the transmission matrix. Within the Random Matrix Theory approach, 
these traces are random variables whose probability distribution depends on the symmetries of the 
system. For the case of broken time-reversal symmetry, we use generalizations of Selberg's integral 
and the theory of symmetric polynomials to present explicit closed expressions for the average value 
and for the variance of T„ for all n. In particular, this provides the charge cumulants {{Qn)) of all 
orders. We also compute the moments {g") of the conductance g — Ti. All the results obtained are 
exact, i.e. they are valid for arbitrary numbers of open channels. 

PACS numbers: 73.23.-b, 05.45.Mt, 73.63.Kv 
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I. INTRODUCTION 

We consider the problem of electronic transport 
through a ballistic quantum dot attached to two ideal 
leads supporting A^i and N2 open channels. Even at zero 
temperature, the electric current as a function of time 
has a fundamentally random nature associated with the 
granularity of charge, and can naturally be character- 
ized by its statisticsi. Its average and variance, for ex- 
ample, are related respectively to the conductance and 
the shot-noise. These quantities have been under exper- 
imental investigation for already quite some time, and 
when the classical dynamics in the dot is chaotic they 
can display characteristic features like weak localization^, 
universal conductance fluctuations'^ and constant Fano 
factor—. More recently, much attention has been devoted, 
both theoretically^ and experimentally^, to higher cumu- 
lants comprising the full counting statistical^, which con- 
tain more refined information about the transport pro- 
cess. 

This problem can be described, in the Landauer- 
Biittiker scattering approach, in terms of the unitary 

r t' 



scattering matrix S 



t r' 



or in terms of a her- 



mitian matrix tt^ ^ xhe transmission matrix t has dimen- 
sion N2 X Ni, where A^i and N2 are the number of open 
quantum channels on the entrance and exit lead, respec- 
tively. The total number of channels is = A^i -I- iV2. 
The dimensionless conductance is given hy g — Tr{tt^), 
and more generally one is interested in the quantities 
= TV[(ttt)"]. 

A possible way to model chaotic cavities is to assume 
S to be drawn at random from an appropriate ensem- 
ble, determined only from the existing symmetries. The 
average value of Tn over such an ensemble is denoted 
(Tn). This is the random matrix theory approach^. It 
treats the eigenvalues of ti^ as correlated random vari- 
ables, and successfully predicts all the universal results 
already mentioned. However, even though experiments 
may be done with relatively small channel numbers, so 
far most explicit expressions are valid only in the asymp- 



totic limit Ni,N2 ^ or to the first few terms in a 
perturbative expansion in 

A few particular results exist that are valid for gen- 
eral Ni,N2, like the average and variance of g^iH and 
the density of transmission eigenvalues for equal leads 
A^i = N2.^' More recently, results were found for: {T2) 
and (Tg^^ii^; some nonlinear statistics like variance of 
shot-noise, skewness and kurtosis of conductance^^; the 
density of eigenvalues for unequal leads fii as well as an 
expression for (Tn). The purpose of this paper is to estab- 
lish general exact results, of which several of the above- 
mentioned ones are particular cases, for systems without 
time-reversal symmetry. It is well known that within this 
symmetry class leading-order perturbative results some- 
times turn out to be exact (e.g. for conductance and for 
the form factor of closed systems). However, this is by no 
means usual, and most transport statistics are in general 
not identical with their N 00 asymptotic limit. 

We shall start by computing (Tn). Our formula is dif- 
ferent from the one in Ref. [isl and computationally supe- 
rior. We also present the variance of Tn and all moments 
of the conductance, {T{^). These advancements are made 
possible by making a connection with the theory of sym- 
metric polynomials and generalizations of Selberg's inte- 
gral. 

We must remark that the quantities more readily ac- 
cessible to measurement are the charge cumulants {{Qn)), 
which quantify the fluctuations in the amount of charge 
transmitted over an interval of time. - They are related 
to {Tn) according to the generating function 



n=l 



i-iy 



■{Tn.){e^-lY 



(1) 



Our exact expression for {Tn) therefore provides cumu- 
lants {{Qn)) of arbitrary order. 

For broken time-rcveral systems, random matrix the- 
ory predicts that the non-zero eigenvalues {Ti,...r/Vj} 
of tt^ behave statistically like random numbers dis- 
tributed between and 1 according to the joint prob- 
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ability distribution^iii 

V{T) =JV\A{T)\^Y[T°', 



(2) 



where A(r) = Y[i<ji'^i ~ ^j) is the Vandermonde deter- 
minant, a = N2 — N1 (we assume iVi < N2) measures the 
asymmetry between the leads and TV is a normalization 
constant. Averages are therefore obtained as 



(/(T)) = / r{T)f{T)dT, 



(3) 



where C is the hypercube [0, . By integrating out all 
but one of the variables one obtains the eigenvalue den- 
sity p(T). An exact expression for p{T) can be written 
down using Jacobi polynomials ) 

Ni-l 

p{T) = T" 5] (a + 2j + l)[P^''-°\l - 2T)]2, (4) 
i=o 

and hence (T„) = Ni T''p{T)dT. Our resuh ^ may 
be seen as an explicit solution to this integral. 



II. GENERALIZATIONS OF THE SELBERG 
INTEGRAL 

We start by fixing some notation (a basic reference 
is dl]). A non- increasing sequence of positive integers 
A = (Ai, A2, . . .) is said to be a partition of n — 
This is indicated by A h n. The number of parts in A is 
called its length and denoted by ^(A). The set of all 
partitions of n is P{n). With every partition A we can 
associate a monomial x'^ = x^^Xn^ • 



■ ■ ■ . The symmetric 
polynomial in k variables m^-^'^ (x) is the sum of all distinct 
monomials obtainable from x^ by permutation of the x's. 
For example, take n — k — 2. Then m(2)(a;) — x\-\- x\ 

and m(i_i)(x) — x\X2- The set of functions {iri\' (x\ A G 
P(n)} forms a basis for the vector space of homogeneous 
symmetric polynomials of degree n in A; variables. 

Another basis for this space is composed by Schur func- 



tions s^^ix). 



If A h n < fc, 



,f{x) 



A(x) 



-det 



(5) 



Let us use again n = fc = 2 as an example. We then 
have S(2)(a;) = x\X2 and S(i — x\^x\^x\X2 - These 
two basis are of course related by a linear transformation, 

*A^^ = S/jhn ■ The matrices K are called Kostka 

matrices. 

Let us from now on take k — Ni, i.e. consider only 
polynomials in A^i variables, and no longer write any su- 
perscript. Let 



[x\x = 



n 



{x + Xi-i)\ 
{x-i)\ 



(6) 



be a generalization of the usual rising factorial 

{x)a=x{x + l)---{x + a-l). (7) 

Our favorite examples are now [a;] (2) = x{x + 1) and 
[2^] (1.1) = — 1). Let us also define the function 



n,<,(A,-A,+j-i)- 



(8) 



This is equal to the product of all hook lengths of the 
partition A,^- but that is not relevant here. Kanckc — 
and Kadell^^ have proved a generalization of Selberg's 
integral which implies that 



..(T)7'(r)dr = I^^ilfI^. (9) 

C [J>l\\n\ 



Other generalizations arc known. For example, HuE^i^ 
(see also [2^) has shown that 



sx[T)s^{T)r{T)dT = {sx){s^){N + (10) 

ic 
where 



[N -I- j + \+ p,){N -i- ]) 



III. RESULTS 



A. Average of Tn 



Making use of ([9]) one can in principle obtain the aver- 
age value of any symmetric polynomial of the transmis- 
sion eigenvalues. For example, it is knownlS, that 



ri-l 



T„ = m(„)(T) = > (-l)''s(„_p,iP)(T), (12) 



p=0 



where the notation V means that the number 1 appears p 
times. Substituting into ([9]) gives the remarkably simple 
result for the counting statistics 



{rn) 



n-1 
p=0 



where 



G 



n-l\ {Ni - p)r,{N2 - p)r 
P 



{N-P)n 



(13) 



(14) 



The sum (jl3p allows in principle the calculation of any 
linear statistic. It has the merit of having only n terms, 
in contrast with the result in [l^ where the number of 
terms in the sum grows with Ni,N2. In Fig.l (top) we 
plot (Tn) /N for a few values of n, as a function of A^i for 
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FIG. 1; Top: some of the counting statistics for 7^2 — 20. The 
value at the maximum iVi « N2 is close to (^„")/4". Bottom: 
Comparison between exact result (I13|l and asymptotic limit 
^ for N2 = 1. 



N2 = 20. All curves have a maximum around Ni ^ N2, 
where they are approximately equal to (^^") /4". 

Let us denote by (7^) 00 the asymptotic form of {Tn) 
when both channel numbers are large, Ni,N2 » 1. 
These functions are given by2. 



{Tn)oo = 



m=0 



-1) 



2m 
m 



1 



(15) 



where we have defined the finite variable ^ = NiN2/N'^- 
Interestingly, a derivation of this expression directly from 
does not seem to be trivial. Notice that (7^)oo = 
(^^')/4" for Ni= N2. In Fig.l (bottom) we compare the 
exact result ([T^ with the asymptotics ^T5\i for N2 — 1, 
in which case (7^) = Ni/{Ni + n). We can see that 
these quantities may differ significantly if both channel 
numbers are close to unity, specially for higher values of 
n. 

We can obtain from (fT5|) the set of charge cumulants 
{{Qn)), which are given by 



{{Qn)) — ^ ^ Cn,m{%n) 1 
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FIG. 2: (color online) Charge cumulants {{Qn}) as functions 
of A^i for iV2 = 3 (top) and 7V2 = 10 (middle). In the bottom 
panel we plot the difference -D„ — {{Qn}) — {{Qn))oa between 
the exact result and its asymptotic expression for equal leads. 



factor (A^2 — Ni)'^ in general. Except for the second one 
(shot-noise), all {{Qn}) oscillate and even change sign 
as functions of iVi for fixed 7V2. The asymptotic form 
of these quantities for large channel numbers, which we 
denote {{Qn))oc, can be derived by using (7^)oo instead 
of {%,) in dill). That yields^ 



{{Qn))c 



(16) 



n-l 
m=0 



. (2m)! 
(m - 1)!' 



(18) 



with the coefficients 

m— 1 



j=0 ^ 



(17) 



Figure 2 shows some plots of {{Qn))- The odd ones are 
always zero for equal leads, A^i = A^2, and contain a 



where Sn,m are the Stirling numbers of the second kind. 
In the last panel of Fig. 2 we compare the exact and 
asymptotic results for A^i = iV2 = A^. In contrast to {Tn), 
for which the asymptotic value is always above the exact 
one, here we see that the difference {{Qn)) — ((Qn))oo 
oscillates as a function of N, with a number of zeros that 
increases with n. 
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FIG. 3: Variances var(7^) as functions of A^i for N2 ~ 20. 

B. Variance of Tn 

We now turn to the calculation of var(7^i). First, we 
note that, analogously to p^ . we can write 

ri-l 

- E (-ir+^S(«-p,l-)(r)s(„-,,l')(r). (19) 

Using PH]) we see that var(T„) = (T^^) - (T^i)^ is given 
by 



var(7;,) 



n-1 

^ (A^ + 2n - w - g - 1)(7V - w - g - 1) ■ 



(20) 

This generalizes to n > 1 the well known result for con- 
ductance fluctuations, var(Ti) = N^N^/N^iN"^ - 1). Of 
course, we could just as easily have computed the value 
of (TnTm) for any n,m. 

In Fig. 3 we plot var(7^) as a function of A^i for the 
case N2 = 20, for a few values of n. The function in- 
creases mildly with n, and is monotonic in Ni. They all 
saturate at a finite value as A^i — > 00. For iVi — N2 the 
many-channels limit can be obtained explicitly by other 
means-- and is given by 

(2n-l)r(n + l/2)r(n-l/2) 

= 8nn{T{nW " ^''^ 

This quantity starts at the value 1/16 (variance of the 
conductance), and saturates at l/47r for large n. 



C. Moments of the conductance 

It is known that for large channel numbers the proba- 
bility distribution of the conductance g = Ti approaches 
a Gaussian, while for small Ni,N2 there are deviations. 
The last quantity we shall compute is the set of all mo- 
ments ((7") characterizing this distribution. Forrester- 
has proved (see also [23,]) that 



2 A 



(22) 




FIG. 4: Top: Even cumulants ({g^")) of the conductance for 
equal leads Ni — N2 = N (odd cumulants are identically 
zero). Middle: {{g")) as a function of iVi for N2 =20. Bot- 
tom: Decay of ((ff")> with n for 7V2 = 2Ni = 20. Notice the 
logarithmic scale in all cases. 



where we have an infinite sum over partitions. Let us 
perform the change of variables x = y/a and take the 
limit a 00. It is well known^^ that, if A h n, then 



[a]A«a", SA(2;/a)«(^) 



IJA 



Hx 



(23) 



Using limQ^oo(l — yT/a)""" = e^^ on the left hand 
side and the simplifications above on the right hand 
side, equation ([^ provides a generating function for the 
quantities we are looking for. The result is 



[iVl]A[iV2]A 



(24) 



From these moments one may compute the corre- 
sponding cumulants ((g")). For a Gaussian distribution. 
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{{g")) = for n > 2. From the symmetry of the distri- 
bution ([2]), aU odd cumulants with n > I vanish iden- 
tically if A^i — N2 (in general, they contain the factor 
(A^i — -^^2)^)- Figure 4 (top) presents some of the even 
cumulants ((g^")) for Ni = N2 = N. We see that ((5^)) 
saturates at the expected value 1/16. Although the oth- 
ers seem to have somewhat singular behaviour for N < n, 
only integer values of N are actually physical. They all 
decrease rapidly after N > n, and already for iV = 6 
the distribution function of the conductance is Gaussian 
to a good approximation. In Fig. 4 (middle) we see the 
first few cumulants as functions of Ni for 7V2 = 20, while 
Fig. 4 (bottom) shows how fast the value of ((g")) decays 
with n, for N2 = 2Ni = 20. 

IV. CONCLUSIONS 

We have extended the random matrix theory of quan- 
tum transport in ballistic chaotic cavities connected to 
two ideal leads, by obtaining exact results for several sta- 
tistical quantities. All our formulas as expressed in sim- 
ple closed form. Even though we consider systems with 
broken time-reversal (TR) symmetry, we have seen that 
exact results are in general different from their large-N 
asymptotic limit. 



Unfortunately, the calculations presented here are not 
immediately generalizable to TR symmetric systems, for 
which we expect the difference between exact and asymp- 
totic results to be more important. This is because the 
role played here by Schur functions is played in that case 
by zonal polynomials, and much less is known about the 
analogues of Kostka matrices. This certainly deserves 
further investigation (actually, relation has been 

proven22 for other symmetry classes, so the moments of 
the conductance can in fact be obtained in those cases). 
Also, it would be interesting to obtain systematic exact 
results for more general nonlinear statistics. 

After this work was completed, another paper 
appeared^l in which the probability distribution of con- 
ductance is studied for broken TR. In particular, they 
obtain an efficient recurrence relation for the cumulants 
of g. It would be important to clarify what are the rela- 
tions between those results and the present work. 
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